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Q \ Abstract 

We analyze the nonlinear relaxation of a complex ecosystem composed of many in- 
teracting species. The ecological system is described by generalized Lotka-Volterra 

k> , equations with a multiplicative noise. The transient dynamics is studied in the 

\^ ' framework of the mean field theory and with random interaction between the 

5r ! species. We focus on the statistical properties of the asymptotic behaviour of the 

time integral of the i-th population and on the distribution of the population and 
of the local field. 



1. INTRODUCTION 

Systems of interacting biological species evolve through a dynamical complex process that 
can be conveniently, described, within relatively short time scales, by generalized Lotka- 
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vestigation, especially in the case of a great number of interacting species. 
Nevertheless some analytical approximation for the mean field interaction between the 
species as well as numerical simulations give some insight into the behaviour of complex 
ecosystems ufl'Q. Basic elements of a Lotka-Volterra model are the growth parameter and 
the interaction parameter. For a large number of interacting species, it is reasonable, as 
a phenomenological approach, to choose these parameters at random from given probabil- 
ity distributions. Within this type of representation, the dynamics of coevolving species 
can be characterized by statistical properties over different realizations of narameter sets. 
Though the generalized Lotka-Volterra model has been explored in detail a, it seems that 
a full characterization, either deterministic or statistical, of the conditions under which a 
population extinguishes or survives in the competition process has not been achieved H'Q. 

In this paper we analyze the role of the noise on the transient dynamics of the ecosystem 
of many interacting species in the presence of an absorbing barrier, i. e. extinction of 
the species. Two type of interaction between the species are considered: (a) mean field 
interaction, and (b) random interaction. We focus on the asymptotic behaviour of the time 
integral of the i-th population and on the distributions of the population and of the local 
field, which is the total interaction of all species on the i-th population. 

By introducing an approximation for the time integral of the average species concentration 
M{t) we obtain analytical results for the transient behaviour and the asymptotic statistical 
properties of the time average of the i-th population. We find that for a very large number of 
interacting species the statistical properties of the time average of the i-th population process 
are determined asymptotically from the statistical properties of the Wiener process. At the 
critical point and around the stability-instability transition, for mean field interaction, the 
system goes from a purely long time tail behaviour (namely M{t) ~ \/t to a new long 
time modified regime M{t) ~ \/te^*. Specifically for random interaction we find that the 
local field and the cavity field 0, which is the total interaction of all species on the i-th 
population when this population is absent, are different in absence of external noise while 
overlap quite well in the presence of the noise. This behaviour, which is very different 
from the analogue spin glasses problem in statistical mechanics, is reminiscent of a phase 
transition phenomenon. It suggest that, becuase all population are positive and can grow 
during the dynamical process of the ecosystem, each population play an important role on 
the total interaction between the species. 

The paper is organized as follows. In the next section we describe the model. The mean 
field and the random interactions are considered in section 3 and 4. Our results are discussed 
in section 5. 



2. THE MODEL 

We consider an N-species generalization of the usual Lotka-Volterra system with a Malthus- 
Verhulst modelization of the self regulation processes for a fully connected ecological net- 
work. Therefore the Ito stochastic differential equation describing the dynamical evolution 
of the ecosystem is 



dni{t) 



^+9) -ni{t)+^Jijnj{t) 



ni{t)dt + ^/eni{t)dwi^ i = l,...,N (1) 
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where nj(i) > is the number of elements of the i-th species. In Eq.(|I|) 7 is the growth 
parameter, the interaction matrix Jij modehzes the interaction between different species 
{i 7^ j) and Wi is the Wiener process whose increment dwi satisfies the usual statistical 
properties 



< dwi{t) > 



0; 



< dwi{t)dwj{t ) > = 5ij6{t — t )dt. 



(2) 



Our ecosystem is composed of A^ = 1000 species. We consider all species equivalent so 
that the characteristic parameters of the ecosystem are independent of the species. The 
random interaction with the environment (climate, disease, etc..) is taken into account by 
introducing a multiplicative noise in the Eq.(|l|). The solution of the dynamical equation 
Eq.(|l|) is given by 



ni{t) 



ni{0)exp 6t + y/ewi{t) + /q dt' Y.j^i Jijnj{t') 



1 + 7ni(0) Jo* dt'exp 5t' + ^ew^{t') + J^*' dt" Y.j^^ Jijnj{t"] 



(3) 



We consider two different types of interaction between the species: (a) a mean field ap- 
proximation with a symbiotic interaction between the species; (b) a random interaction 
between the species with different types of mutual interactions: competitive, symbiotic and 
prey-predator relationship. 



3. MEAN FIELD APPROXIMATION 

We consider a mean field symbiotic interaction between the species. As a consequence the 
growth parameter is proportional to the average species concentration 



and the stochastic differential equation Eq.(|lj) becomes 



(4) 



drij 



Jm + 7 + - ) n. 



n,- 



dt + y/eriidwi 



(5) 



In the limit of a large number of interacting species the stochastic evolution of the system 
is given by the following integral equation 



where 



M(t) = -Y^\n{l + n,(0) /■*dt'e^*^(*')+^*'+v^'"'(*')] , 

^W = ^E fdt'mit') = fdt'm{t'). 
1\ , Jo Jo 



(6) 



(7) 



is the time integral of the site population concentration average. We introduce an approx- 
imation of Eq.(g) which greatly simplifies the noise affected evolution of the system and 
allows us to obtain analytical results for the population dynamics. We note that in this ap- 
proximation the noise influence is taken into account in a nonperturbative way, and that the 
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the statistical properties of the process Wmax{t) = supo<f/<t z«(t')) where w is the Wiener 
process. Starting from the following approximated integral equation for M{t) 



M{t) ~ — ^ In ('l + ni(0)ev^"'™'"-« f 



^^,^JM(t')+^t' 



(8) 



it is possible to analyze the role of the noise on the stability-instability transition in three 
different regimes of the nonlinear relaxation of the system: (i) towards the equilibrium 
population (7 > 0), (ii) towards the absorbing barrier (7 < 0), (iii) at the critical point 
(7 = 0). Specifically at the critical point we obtain for the time average process M{t)/t as 
a dominant asymptotic behaviour in the stability region (namely when J < 1) 



M{t) 



1 



t \l-J 

and in the instability region (namely when J > 1) 



2e 1 



(9) 



M{t) 



={lnK(0))> 



2-Kexp 



Vt 



Vt 



(10) 



We obtain for the case (i) an explicit expression of the transition time tc as a function of 
the noise intensity (e), the initial population distribution (ni(0)) and the parameters of the 
system (7, J) 




2-irj 



+ In 1 + 



7 



(J-1) 



-{HMO))) 



1/2 



27r7 



:ii) 



For the cases (ii) and (iii) we obtain two implicit expressions in terms of exponential and 
error functions of the same quantities: tc, e and the system paramteraJ. The transition time 
increases from 7 > to 7 < according to the following inequality 



{tc)'y<0 > {tc)^=0 > (ic)7>0- 



(12) 



This means that when the interaction between the species prevails over the resources, the 
presence of a hostile environment (7 < 0) causes a late start of the divergence of some 
population (i.e. the instability). The noise forces the system to sample more of the available 
range in the parameter space and therefore moves the system towards the instability. The 
effect of the noise is to make unstable the system earlier than in the deterministic case 
(e = 0). If we raise the intensity of the noise and keep fixed the initial distribution, we 
obtain the same effect of the enhancement of the variance of the Gaussian initial distribution 
of the population for moderate values of noise intensity (namely e = 0.1). For high values 
of noise intensity (namely e = 1) we strongly perturb the population dynamics and because 
of the presence of an absorbing barrier we obtain quickly the extinction of the populations. 



4. RANDOM INTERACTION 

The interaction between the species is assumed to be random and it is described by a 
random interaction matrix Jij, whose elements are independently distributed according to 



P{J^ 



^J) 



27raj 



-.exp 



where J is the interaction strength and 
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-1 



N' 



(13) 



< Jij >-- 



0, 



\ JijJji ^- 



0. 



(14) 



With this asymmetric interaction matrix our ecosystem contains 50% of prey-predator in- 
teractions (namely Jij < and Jji > 0), 25% competitive interactions {Jij < and Jji < 0) 
and 25% symbiotic interactions {Jij > and Jji > 0). The initial values of the populations 
nj(0) have also Gaussian distribution 



P(n) 



1 



V2^ 



exp 



{n— < n >)^ 



a^ = 0.01, and < n >= 1. 



(15) 



The strength of interaction between the species J determines two differnt dynamical be- 
haviours of the ecosystem. Above a critical value Jc the system is unstable, this means that 
at least one of the populations diverges. Below the critical interaction strength, the system 
is stable and reaches asymptotically an equilibrium state. For our ecosystem this critical 
value is approximately J = 1.1. The equilibrium values of the populations depend both on 
their initial values and on the interaction matrix. If we consider a quenched random inter- 
action matrix, the ecosystem has a great number of equilibrium configurations, each one 
with its attraction basin. For vanishing noise (e = 0), the steady state solutions of Eq.(||) 
are obtained by the fixed-point equation 



(7 - nj + hi)ni = 



(16) 



where 



hi = ^Jijnj{t) 



(17) 



is the local field. For a great number of interacting species we can assume that the local 
field hi is Gaussian with zero mean and variance a^. =< hf >= J^ < ^f > 



P{h) 



1 



2^^i 



exp 



h} 



2-L 



(18) 



The solutions of Eq. (|lq) are 



and 



Uj = 0, i. e. extinction 



"-j = (7 + hi)Q{l + hi), Ui > 0, 



(19) 



(20) 



where Q is the Heaviside unit step function. From this equation and applying the self 
consistent condition we can calculate the steady state average population and its variance. 
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< rii > 



{{l + K)Q{^ + hi)) 



^^< 



a^^exp 



7 
2-1 



+ 



7 



^<^ 



1 + erf 




(21) 



and 



< nf > 



1 + erf 



'l' + < 




^7 /2< 



7 



2-1 



(22) 



For an interaction strength J = 1 and an intrinsec growth parameter 7 = 1 we obtain: 
< Ui >= 1.4387, < nf = 4.514, and o"^^ = 2.44. These values are in good agreement with 
that obtained from numerical simulation of Eq. (||). The choice of this particular value 
for the interaction strength, based on a preliminar investigation on the stability-instability 
transition of the ecosystem, ensures us that the ecosystem is stable. 

The stationary probability distribution of the populations is the sum of a delta function 
and a truncated Gaussian 



exp 



Pirii 



neiS{ni) + @{ni)- 






2^,Pal. 



(23) 



In Fig.(l) we report the stationary probability distribution of the population densities, 
without the extinct species, in comparison with the computer simulations for systems with 
N = 1000 species and for an interaction strength J = 1, and 7 = 1. 

As in the previous case we focus on the statistical properties of the time integral of the 
i-th population Ni{t) 



(24) 



N,{t) = / dt'riiit'), 
Jo 

in the asymptotic regime. Prom Eq. (pi) we have 



Ni{t)= In 



1 + nj(0) / dt'exp 






(25) 



In Eq. (^) the term ^ ■ JijNj gives the influence of other species on the differential growth 
rate of the time integral of the i-th population and represents a local field acting on the ith 
population U'LI 



hi =^JijNj{t) = Jrji. 



(26) 



We use the same approximation of the Eq.(|^ and, after differentiating, we get the asymp- 
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Ni{t) ~ In 



n. 



jQ\^y^U)^ax^it) + JVmaXi(.t) f ^^'gl*' 



(27) 



where Wmaxiit) = supo<t'<tw{t') and r]maxi{t) = supo<,t'<tri{t')- The Eq.(|^ is vahd for 
7 > 0, that is when the system relaxes towards an equihbrium population and at the 
critical point. Evaluating Eq. (|27|) for 7 > 0, after making the ensemble average, we obtain 
for the time average of the i-th population iVj 



{N,) ~ - [iV^Vet + Int + (In [n,{o)]) 



7 = 0, 



and 



{Ni)^ 



NWet + (7 + A^^ + ( In 



7 



7>0, 



(28) 



(29) 



where A'^ and Njj are variables with a semi-Gaussian distribution □ and N^j must be deter- 



mined self-consistently from the Eq. (|2q). These asymptotic behaviours are consistent with 
those obtained using a mean field approximation. We obtain in fact the typical long time 




Figure 1: The stationary probability distribution (Eq.(p3|)), without the extinct species, 
in comparison with the histograms arising from numerical simulations (open circles). The 
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tail behaviour (t~^") dependence, which characterize nonhnear relaxation regimes when 
7 > 0. Besides the numerical results confirm these analytical asymptotic behaviours of iVj 
□. When the system relaxes towards the absorbing barrier (7 < 0) we get from Eq. (|2~ 
the long time regime 



m 



{N^) ^ I 



ln(ni(0))+ln 



^£^lt'+Ww,{t')+Jm{t') 



(30) 



In this case the time average of the i-th population {Ni) is a functional of the local field 
and the Wiener process, and it depends on the history of these two stochastic processes. 
We have also analyzed the dynamics of the ecosystem when one species is absent. Specifically 
we considered the cavity field, which is the field acting on the i-th population when this 
population is absent 0. In Fig. (2) we report the probability distributions of the local and 
of the cavity fields obtained by our simulations after a time t = 100 (expressed in arbitrary 
units) in absence of external noise and for two species (namely species 1 and 33). We 
note that the probability distributions of the cavity fields differ substantially from that of 
the local fields for the same species unlike the spin glasses dynamics, where the two fields 
coincide. We calculate also the same quantities in the presence of the external noise. The 
results of our simulations are reported in Fig.s (3) and (4). The effect of the external noise 
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Figure 2: The probability distribution of the cavity fields P{hc) (open circles) and of the 
local fields P{hi) (black circles) for the species 1 and 33 after time t = 100. The system 
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is to overlap the two fields in such a way that for some particular species they coincide. 
Specifically this happens for the species 1 (see Fig. (3)). For the species 33 we obtain a 
partial overlap (see Fig. (4)). 

We found this interesting phenomenon, which is reminiscent of the phase transition phe- 
nomenon, for some populations. The main reasons for this behaviour are: (i) all the popu- 
lations are positive; (ii) the particular structure of the attraction basins of our ecosystem; 
(iii) the initial conditions, which differ for the value of one population, belong to diffeent 
attraction basins. Some populations have a dynamical behaviour such that after a long time 
they influence in a significant way the dynamics of other species. While in the presence 
of noise all the populations seem to be equivalent from the dynamical point of view. We 
found also that for strong noise intensity (namely e = 1) all species extinguish on a long 
time scale (t ~ 10^ a. u.). Whether extinction occurs for any value of noise intensity or not 
is still an open question, because of time-consuming numerical calculations. 



5. CONCLUSIONS 

0.2 




Figure 3: The probability distribution P{hi) of the local (black circles) and of the cavity 
(open circles) fields for the species 1 after time t = 100, in the presence of external noise. 
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We studied a stochastic model of an ecosystem of N interacting species. By means of an 
approximation of the integral equation, which gives the stochastic evolution of the system, 
we obtain analytical results reproducing very well almost all the transient. We investigate 
the role of the noise on the stability-instability transition and on the transient dynamics. For 
random interaction we obtain asymptotic behaviour for three different nonlinear relaxation 
regimes. We obtain the stationary probability distribution of the population, which is the 
sum of two contributions: (i) a delta function around n = for the extinct species and 
(ii) a truncated Gaussian for the alive species. When we switch on the external noise 
an interesting phenomenon is observed: the local and the cavity fields, whose probability 
distributions are different in the absence of noise, coincide for some populations. This 
phenomenon can be ascribed to the peculiarity of the attraction basins of our ecosystem. 
We have also investigated the overlap between the asymptotic values of the populations and 
the eigenvector of the interaction matrix with the maximum eigenvalue and we have not 
found any ordering regime phenomenon like in the spin glasses system. A more detailed 
investigation concerning the probability distribution of the populations and the local fields 
in the presence of noise is the subject of work in progress. 




Figure 4: The probability distribution P(/i33) of the local (black circles) and of the cavity 
(open circles) fields for the species 33 after time t = 100 in the presence of external noise. 
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